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ABSTRACT 

Recent estimates of magnetic field strengths in T Tauri stars yield values B = 1-4 kG. 
In this paper, I present an upper limit to the photospheric values of B by computing 
the equipartition values for different surface gravities and effective temperatures. The 
values of B derived from the observations exceed this limit, and I examine the possible 
causes for this discrepancy. 
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1. Introduction 

Magnetic fields are believed to play a fundamental 
role in the structure and evolution of T Tauri stars. 
Unfortunately, their detection and measurement in 
these objects is a very difficult task. 

Early attempts to detect magnetic fields in T Tauri 
stars through Zceman polarization were unsuccessful 
(Johnstone & Penston 1986, 1987), and only recently 
such detections are forthcoming through the measure- 
ment of Zeeman broadening of photosphcric lines in 
a handful of objects ( Basri et al. 1992; Guenther 
1997). These measurements yield disk-averaged, pho- 
tospheric field strengths of order a few kG. 

These observations do not provide a direct mea- 
surement of B\ instead, what is measured is the dif- 
ference in broadening of magnetically-sensitive and 
insensitive lines (Robinson 1980; Marcy 1982; Saar 
1988). To find the value of B, one has to fit the 
data with a stellar atmosphere model and vary B and 
/ — the fraction of stellar surface that is magnetized — 
until a good match is found. The result of this match- 
ing procedure is the product fB, and independent 
determination of B and / is much more problem- 
atic and sensitive to modeling assumptions (see, e.g., 
Solanki 1992 and references therein). Furthermore, 
when equivalent widths are used in lieu of measur- 
ing the line broadening, the deduced fB is also very 
sensitive to the input model atmosphere (Basri et al. 
1992; Basri & Marcy 1994). 

The purpose of this paper is to present a physi- 
cal constraint on the maximum value of fB; this up- 
per limit on the magnetic field strength is important 
to estimate B from the fB-values obtained from the 
data, and may help to uncover potential errors in the 
modeling assumptions. 

The basic physics is rather well known, and has 
been successfuly applied to the study of the Sun. 
Magnetic fields are excluded from closed circulation 
patterns when the field lines are perpendicular to the 
axis of rotation of the fluid (see, e.g., Parker 1979, 
ch. 16 and references therein); as a result of this 
general principle the solar convection sweeps the field 
into the downdraft regions of the eddies. Convective 
heat transport is strongly inhibited within the field 
(Biermann 1941), and, through the so-called supera- 
diabatic effect (Parker 1978), the downflowing gas 
further cools and evacuates the magnetized g 
a result, the external gas pressure is larger than that 
inside the tube, and the flux tube shrinks until the to- 



tal internal pressure (thermal plus magnetic) balances 
the external pressure. Thus, the solar magnetic field 
is broken up and compressed into highly evacuated 
flux tubes. 

Although we cannot resolve features on the surface 
of other stars, it is a reasonable inference that, in com- 
plete analogy with the Sun, the magnetic field in stars 
with outer convection zones — or convective through- 
out, as in the case of T Tauri stars — is concentrated 
into discrete flux tubes, which are in pressure equilib- 
rium with the unmagnetized gas. If this is true, then 
an upper limit on B can be derived if the pressure of 
the unmagnetized gas is known and the gas pressure 
inside the flux tubes is ignored 1 . These arguments 
have been extensively applied to active dwarfs (see, 
e.g., Saar 1996), but, so far, they have never been 
used to study the magnetic fields in T Tauri stars. 

In §2 I derive an expression for the maximum 
cquipartition value of B as a function of the stellar 
and atmospheric parameters, and in §3 I compare the 
results from §2 with current observational estimates. 
My conclusions follow in §4. 

2. Equipartition Fields 

Consider a straight flux tube that is initially per- 
pendicular to the photospheric surface at some fidu- 
cial height. The condition of lateral pressure equilib- 
rium, i.e., the flux tube neither expands nor contracts, 
is given by 



where AP is the difference between the external and 
internal non-magnetic pressures, Pi and P e , respec- 
tively, and the external medium is unmagnetized. 
At any height, an upper limit to B is given by the 
equipartition field oc (Pe) 1 / 2 , and because tempera- 
ture and pressure decrease with height in the photo- 
sphere, the maximum detectable field is the equipar- 



1 Note that this argument does not apply to stars with outer ra- 
diative zones, where there are no known mechanisms to isolate 
and concentrate the field. The Ap stars are a prime example, 
with measured fields as high as 30 kG in the case of HD 215441 
(Babcock 1960), a value much larger than that obtained by as- 
suming pressure balance. Therefore, these fields arc likely to 
thread the entire surface, and probably are fossil, rather than 
dynamo-generated, fields. In the case of T Tauri stars, after 
the onset of convection the fossil field is quickly destroyed by 
the resulting turbulent difussivity. Kiiker & Riidiger (1997) 
estimate that in a typical TTS the fossil field is destroyed in 
~ 300 yr after the onset of convection. Therefore, the magnetic 
field in TTS must be dynamo-generated. 
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tition field at the optical depth where the continuum 
is formed (r = 2/3), 



B cq = [8^P e (r-2/3)] 1/2 . 



(2) 



In principle, P e includes the contributions from 
both thermal pressure and ram pressure due to con- 
vective motions (Pt and P v , respectively) because in a 
convective atmosphere the magnetic field outside dark 
spots is concentrated where the convective downdrafts 
are located (see, e.g., Stein ct al., 1992). However, it 
is easy to show that Pt 3> P v , and the argument is 
the following. 

The ratio P v /Pt can be written as 



Pt 



1 



■7 



(3) 



where 7 is the ratio of the specific heat capacities, v is 
the average speed of a (turbulent) convective element, 
and v s is the adiabatic speed of sound. Although re- 
cent simulations of compressible convection (Catta- 
nco & Malagoli, 1992) show that horizontal surface 
flows can be intermittently transonic, the downflows 
are subsonic, and, because 7 < 5/3, it follows that 
Pv/Pt < 5/9. A more stringent limit on this ratio is 
given by estimating v /v s using the mixing- length the- 
ory of convection, which gives (see, e.g., Cox & Giuli 
1968, eq. 14.64) 



v 

— = 0.4 , 
v s Vl-4L Q 



1/3 



T, 



V 4000K 



-1/2 



(4) 

where L*, T c g, and p are, respectively, the stellar 
luminosity, effective temperature, and photospheric 
mean molecular weight. Therefore, for the stellar pa- 
rameters of a typical TTS, P v /Pt ~ 0.2, and by ne- 
glecting P v the error introduced in B eq is ~ 10% — 
much smaller than the observational errors (see be- 
low). Henceforth I will neglect the contribution of 
ram pressure 2 to _B oq , and use P e = Pt- 

The value of P e (r = 2/3) follows from the condi- 
tion of hydrostatic equilibrium, and to a very good 
degree of approximation (generally, better than a fac- 
tor of 2 or 3; see, e.g., Cox & Giuli 1968), one can 
write 

2 9 



P e (r = 2/3) = -J, 

6 K 



(5) 



2 The same result, essentialy, applies to the Sun, where the 
Maxwell stresses at the surface are an order of magnitude larger 
than the Reynolds stresses. 



where k is the Rosseland mean opacity. The gas tem- 
perature T at t = 2/3 is assumed to be the effective 
temperature of the star, T e g, and, because k is a func- 
tion of density p and T, eq. (5) is an implicit equation 
for p(T = T off ) for a given T cff , g, and n(p, T eff ). Thus, 
one finally obtains 



eq 



167T g 
~ K 



1/2 



(6) 



I have computed B eq for different values of T e g and 
g by solving eq. (5) using the Alexander & Ferguson 
(1994) Rosseland mean opacities, and using eq. (6). 
To gauge the magnitudes of the errors introduced by 
the approximation used to derive eq. (5), I have also 
computed B cq by obtaining P e (T — T c g) from the de- 
tailed stellar- atmosphere models by Guenther et. al 
(1992) and Allard & Hauschildt (1995). The com- 
parison between the two sets of B cq so obtained is 
presented in Table 1. 

The results in Table 1 show that for spectral type 
K5, or earlier, the approximation used in eq. (5) re- 
sults in an underestimate of the field by less than 
10%. On the other hand, for later spectral types the 
field is overestimated by, at most, a factor of 2, but 
this error decreases with decreasing gravity. The rea- 
son for this larger discrepancy at cooler temperatures 
is the increase in opacity due to H 2 once the gas 
temperature drops below 3500 K (Alexander & Fer- 
guson 1994), which results in a strong non-monotonic 
behavior of n with optical depth in the outer layers. 

Therefore, for the typical gravities of T Tauri stars 
(log g <J 4.0) and the spectral types currently acces- 
sible to Zeeman-broadcning measurements (spectral 
type K7 or earlier), the equipartition fields predicted 
from eq. (6) are too small by, at most, 10% when 
compared with detailed atmospheric models.. These 
errors are much smaller than the error bars for the 
measured fB's (see Figure 1 below). 

It is important to keep in mind that the measured 
field depends on the geometrical depth where the 
emission originates. Because flux tubes are strongly 
evacuated, the optical depth inside a magnetized re- 
gion is smaller than in the surrounding atmosphere. 
If the flux tube is slender enough, it will be mostly 
transparent to the radiation crossing its walls, and, 
because of its lower density, it may appear as a bright 
feature against the continuum from its surroundings, 
as the network fields in the Sun do. In the solar case, 
the enhanced temperature compensates for the lower 
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Table 1 

Comparison Between B cq from Equation (6) and Detailed Stellar 

Atmospheres 



Sp. Type 


logs 


B eq (kG) a 


A b 


Atmospheric Model 


Sun 


4.4 


1.2 


0.08 


Guenther et al. 1992 


K5 


4.5 


1.8 


0.06 


Allard & Hauschildt 1995 


Ml 


4.0 


1.5 


-0.13 


Allard & Hauschildt 1995 


M5 


4.0 


1.5 


-0.36 


Allard & Hauschildt 1995 


M5 


5.0 


4.8 


-0.52 


Allard & Hauschildt 1995 



a Computed from eq. (6). 

b Fractional difference; the equipartition field computed from a detailed atmospheric 
model is given by B cq [1 + A] . 



optical depth, and the values of B measured in these 
bright points compare well with B cq (T — T c ff)- If the 
Sun is any guidance, then B cq (T = T e g) should be 
also a good approximation to the fields outside spots 
in other stars. 

On the other hand, larger features like sunspots are 
much more evacuated and cooler than slender flux 
tubes, and at a fixed geometrical height the optical 
depth is much smaller inside the spot. Thus, one can 
see deeper into a sunspot, and they indeed appear 
as dimples on the Sun's surface when seen near the 
limb — a phenomenon known as the Wilson depres- 
sion. Because the continuum from a spot originates at 
a larger geometrical depth, the external pressure there 
is larger, and, therefore, the fields measured inside 
sunspots are larger than those from slender features 
like network fields. In the Sun, the Wilson depression 
can amount to a few scale heights, and thus the mea- 
sured fields in sunspots are ~ 2B cq (T — T c g), i.e., 
of order 3 kG (although, sometimes, values as high as 
5 kG have been measured) . In complete analogy with 
the Sun, one should not expect B cq (T = T c g) to be a 
useful limit — within a factor of a few — to the fields in 
stellar dark spots. The equipartition argument still 
applies; however, to derive a useful limit on B it is 
necessary to know the detailed thermal structure of 
the spot, a problem that still lacks a full solution even 
for the Sun. 



3. Comparison with Current Measurements 

At present, only two set of magnetic- field measure- 
ments for T Tauri stars are available 3 : the measure- 
ments by Basri et al. (1992) and Guenther (1997). 
These observations measure the broadening of Fel 
lines in the wavelength range 5000-7000 A; at these 
wavelengths, a dark spot with temperature T e g — AT 
that covers a fraction f s of the surface of a star 
with T c ff contributes a fraction f s B\(T e g — AT)/(1 — 
f s )B\(T c ff) of the total emission, where B\ is the 
Planck function. Using the typical values f s = 0.3, 
AT = 1000 K, and T cfr = 4000 K (see, e.g., Bouvier et 
al. 1993), one finds that starspots contribute less than 
7% to the light in the range 5000-7000 A. Therefore, 
current measurements are insensitive to the magnetic 
fields in the center of dark spots, where B can be 
larger than B eq at T = T e ff. 

Figure 1 is a comparison between the lower limits 
to B derived from the observations (because / < 1) 
and B eq calculated from (6) as a function of T e g and 
different surface gravities. Also shown is the value of 
-B e q for the Sun; comparison of this value with the one 
derived from (6) for a G2 V star gives a measure of 

3 Donati ct al. (1997) detected circular polarization in V410 Tau 
and HDE 283572 (V987 Tau), but they did not derive a value 
of B from their measurements. Guenther & Emerson (1996) 
attempted to measure the fields in Tap 35 and V410 Tau to 
demonstrate the potential of infrared lines in Zccman studies, 
but they obtained upper limits to B that are consistent with 
the measurements of Basri et al. (1992). 
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the typical errors (Table 1) in the constant-g loci in 
Figure 1 for the spectral types currently accessible to 
Zccman-broadening measurements (spectral type K7 
or earlier). Note that these errors are much smaller 
than the observational error bars. 
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K5 MO M5 
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log T eff (K) 



Fig. 1. — The equipartition magnetic field, B cq , as 
a function of effective temperature T e g and different 
values of surface gravity g. Loci of constant g are 
shown as grey lines and are labeled by the value of 
log g. The field strengths derived by Guenther (1997) 
are shown as filled circles (from left to right): T Tau, 
Lk Ca 15, and Lk Ca 16; the values of B from the 
observations by Basri et al. are shown as white cir- 
cles (from left to right): TAP 10 and TAP 35. The 
arrows indicate upper limits. Also shown is the locus 
of i?eq(T e ff) for the main-sequence (dotted line) and 
the value of B eq for the Sun from the Guenther et al. 
(1992) model is indicated by 0. 

There is a troubling trend in Figure 1, which is 
most obvious for the Guenther (1997) data. The field- 
strengths derived by this author are above the main- 
sequence values of B eq . In other words, the only way 
for these objects to have magnetic fields as strong as 
implied by the data is for them to be more compact 
(log g ;> 4.5) than a main-sequence star of the same 
spectral type. Obviously, this is impossible; moreover, 
log g <J 3.5 for classical TTS such as those in the 
Guenther (1997) dataset, and, also, the absorption 
lines form at r < 2/3. 



Furthermore, since I have assumed / = 1 in deriv- 
ing a value of B from the measured fB, these values, 
again, are lower limits to the field implied by the ob- 
servations. Therefore, the measurements of Guenther 
(1997) seem to overestimate B by at least a factor 
of > 2 (2^) 1/2 in the case of Lk Ca 15 and a factor 
> 2.5 (^) 1/2 for T Tau, respectively. 

The measurements of Basri et al. (1992) for TAP 
35 seem compatible with the expected B cq , in particu- 
lar because they assumed \ogg — 4.0 in their models. 
However, if / < 0.65 their measurements would im- 
ply a value of B larger than the equipartition field 
in a main-sequence star of the same spectral type. 
Moreover, they found that if log g = 3.5, instead of 
log g = 4.0, then a stronger field is required to match 
the data; this is inconsistent with the dependence of 
B eq on g at fixed T e g (Figure 1). Therefore, there 
are also problems with the Basri et al. (1992) mea- 
surements, and the implied fields are larger than B eq . 

Another perspective on this inconsistency is pre- 
sented in Figure 2. This figure is a plot of B eq as a 
function of T e gf along the pre-main-squence isochrones 
of DAntona & Mazzitelli (1994); also shown are the 
paths in the (T e g, -B eq )-plane traced by two stars with 
M = 1.5 M Q and M = 0.4 M , respectively, as they 
evolve towards the zero-age main-sequence. These 
results show that, for a given value of the stellar 
mass, B eq increases with age; therefore, if one takes 
« 10 Myr as the age limit for the T Tauri phase, the 
results in Figure 2 imply that in TTS of spectral type 
K7 or earlier B £ 1.5 kG outside dark spots. 

4. Conclusions 

The equipartition field B cq at any level in a stel- 
lar atmosphere is an upper limit to the magnetic field 
strength in a flux tube in pressure equilibrium with 
the surrounding, non-magnetized gas. I have com- 
puted B eq for the physical regime appropriate to T 
Tauri stars, and shown that the field strengths derived 
from current Zeemann-broadening measurements ex- 
ceed this upper limit by large factors, and therefore 
are unphysical. 

What are the causes for this discrepancy? 

Because the derivation of B from the measured line 
widths requires detailed modeling of the emergent line 
profiles, the fact that current models yield B > B eq 
indicates that some piece of physics is missing. To be 
fair to the Zccman-broadening technique, when the 
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Zceman components are resolved this method can be 
very precise (see, e.g., the recent measurements of B 
in e Eri by Valenti et al., 1995), whereas Basri et al. 
(1992) based their analysis on the equivalent widths 
of the lines (due to the faintness of TTS and the high 
signal-to-noise required by this method), and Guen- 
ther (1997) used an autocorrelation analysis without 
benefit of a detailed atmospheric model. However, it 
is hard to imagine how the use of equivalent widths 
rather than line profiles can explain the finding by 
Basri et al. (1992) that a lower surface gravity re- 
quires a larger magnetic field to match the data, when 
the field is supposed to decrease with decreasing grav- 
ity (see eq. [6]). Therefore, the source of the discrep- 
ancy is, most likely, the input physics of the models. 

In particular, current models for TTS use the 
same atmospheric structure for the magnetic and non- 
magnetic parts of the photosphere, whereas at any 
given height the gas pressure in a flux tube has to 
be lower than that of its surroundings for pressure 
equilibrium to obtain. Because the measured field de- 
pends on the geometrical height where the lines are 
formed, and the flux tubes must be highly evacuated, 
using the same atmospheric structure for the magne- 
tized and quiet regions of the photosphere is not cor- 
rect, and is the most likely origin of the discrepancy. 
Therefore, this and other refinements (see Landolfi et 
al. 1989) must be incorporated to derive a more accu- 
rate estimate of B in TTS. In addition, Valenti et al. 
(1995) showed that the Zeeman-broadcning method 
is much more immune to the details of the model 
atmosphere when infrared lines are used. Infrared 
measurements should be pursued also, and one hopes 
that the discrepancy with the results presented here 
will disappear. 
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Fig. 2. — The equipartition magnetic field, B cq , as a 
function of effective temperature T c g along the pre- 
main-sequence isochrones of D'Antona & Mazzitclli 
(1994). The loci of constant age t are shown as dotted 
lines and are labeled by the value of t. The solid lines 
indicate the paths in the (T e g, -B eq )-plane for a 1.5 M Q 
(left) and a 0.4 M Q star (right) as they evolve towards 
the zero-age main-sequence. The other symbols have 
the same meaning as in Figure 1. 
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